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Abstract 

We are interested in the inverse problem of recovering a Robin coefficient defined on some non accessible 
^ part of the boundary from available data on another part of the boundary in the nonstationary Stokes system. We 
[ prove a Lipschitz stability estimate under the a priori assumption that the Robin coefficient lives in some compact 
and convex subset of a finite dimensional vectorial subspace of the set of continuous functions. To do so, we use a 
theorem proved by L. Bourgeois which establishes Lipschitz stability estimates for a class of inverse problems in an 
abstract framework. 
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p ^ Resume Estimation de stabilite Lipschitzienne pour le systeme de Stokes avec des conditions aux limites de 

<' ' ' 



Nous nous interessons a l'idcntihcation d'un coefficient de Robin defini sur une partie non accessible du 
bord a partir de mesures disponibles sur une autre partie du bord dans le systeme de Stokes non stationnaire. Nous 
prouvons une estimation de stabilite Lipschitzienne sous l'hypothese a priori que le coefficient de Robin est defini 
Ch ■ dans un sous-ensemble compact et convexe d'un sous-espace vectoriel de dimension finie de l'espace des fonctions 
continues. Pour ce faire, nous utilisons un theoreme prouve par L. Bourgeois permettant d'etablir des inegalites de 
— i ■ stabilite Lipschitzienne pour une classe de problemes inverses dans un cadre abstrait. 
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' Soit T > 0, C R d , avec d £ N*, un ouvert borne Lipschitzien et connexe tel que <9fi = L; UL UL otlt avec T out = M L, ; 
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et v est la normale exterieure a f2. On considere le systeme de Stokes suivant : 

d t u - Au + Vp = 0, dans (0, T) x SI, 

divu = 0, dans(0,T)xft, 

u =0, dans (0,T) x T u 

d v u-pv = g, sur(0,T)xL , 

d v u-pv + qu = 0, sur (0,T) x Y out , 

u(0, .) = uo, dans £1. 

Le probleme inverse qui nous intcrcssc est le suivant : on cherche a identifier le coefficient de Robin q defini 
sur la partie non accessible du bord T out a partir de mesures disponibles sur Tq pour (u,p) solution du systeme 
([T]). Ce type de probleme inverse apparait naturellement dans la modelisation d'ecoulements biologiques, comme par 



exemple l'ecoulement sanguin dans le systeme cardiovasculaire (voir et ou encore l'ecoulement de Pair dans 
les poumons (voir Q). Nous renvoyons a 0] pour une introduction a la modelisation de l'ecoulement de l'air dans les 
poumons et aux differentes conditions aux limites qui peuvent etre prescrites. La resolution de ce type de probleme 
inverse dans le cas stationnaire a deja ete etudiee dans Q, 0] et Q. Dans les deux premiers travaux, une inegalite de 
stabilite logarithmique est obtenue alors qu'une inegalite de stabilite Lipschitzienne est etablie dans @ sous l'hypothese 
a priori que le coefficient de Robin est constant par morceaux sur T out . Dans chacun de ces papiers, les mesures 
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dp 

intervenant dans les inegalites de stabilite sont la vitesse u, la pression p et la derivee normale de la pression — — sur 

r-. on 
r C To. Le cas du systeme de Stokes non stationnaire a ete aborde dans [5| dans le cas particulier ou le coefficient 

de Robin ne depend pas du temps. L'idee, introduite dans 0] dans le cas de l'equation de Laplace, consiste a etendre 

l'inegalite de stabilite valable pour le probleme stationnaire au probleme non stationnaire en utilisant une inegalitc 

provcnant de la thcoric dcs scmigroupcs analytiques. Cela conduit a faire des mesures en temps infini. 

L'originalite de l'inegalite de stabilite Lipschitzienne presentee dans cette Note est multiple : d'une part, 
nous obtenons une inegalite de stabilite valable pour le systeme de Stokes non stationnaire en temps fini avec un 
coefficient de Robin dependant du temps et d'autre part, l'unique mesure intervenant dans l'inegalite de stablite est la 
vitesse u sur (0,T) x L, avec L C r . De plus, l'ensemble des coefficients de Robin pour lequel l'inegalite de stabilite 
Lipschitzienne est valide est un peu plus general que dans Q : les coefficients de Robin ne sont plus necessairement 
constants par morceaux mais appartiennent a un sous-ensemble compact et convexe d'un sous-espace vectoriel de 
dimension finie de l'ensemble des fonctions continues. Enfin, nous avons besoin d'hypotheses de regularity moins 
fortes sur le bord du domaine et sur le flux g. 

Afin d'etre plus precis, nous introduisons quelques notations. 
Notation 1.1. On note 

L5f((0,T) x T out ) = {qe L°°((0,T) x T out ); 3m >0, q > m p. p. sur (0,T) x r out }, 



et 



C°(0,T;C a pc (T out )) = {q : (0,T) x T out -> R; q\ {0 ,T)xr t £ C°((0,T) x L 2 ) pour 1 < i < N}. 



Le resultat principal de cette Note est resume dans le theoreme suivant : 

Theoreme 1.2. Soit M 6 N*. On considere Vm un sous-espace vectoriel de C°(0,T;Cp C (T out )) engendre par M 
fonctions lineairement independantes et Km un sous-espace convexe et compact de Vm^LT((0, T) x T out ). Soit V C Tq 
une partie ouverte non vide du bord du domaine, Uq € H (CI) tel que div uq — dans ft et g G H 1 ^, T; L 2 (Tq)) tel que 
g(t) est non identiquement zero pour tout t £ (0,T). Soit (iik,Pk) la, solution faible du systeme ([1} avec q = qk £ Km 
pour k = 1,2. Alors, il existe une constante C > telle que 

1 1 9i - 92|U~((0,T)xr o „ t ) < C\\ui - M2||z, 2 ((o,T)xr)- 

La preuve du Theoreme 11.21 est basee sur un theoreme abstrait prouve par L. Bourgeois dans [6( que nous 
rappelons dans le Theoreme 13.41 et repose sur le fait que l'application 

T : L™((0,T)xT out ) -> L 2 ((0,T)xr) 
q -4 u\ r 

ou (u,p) est solution du systeme (TTJ) et avec L C r , est injective, de classe C 1 et sa derivee est egalement injective. 
Notons que le resultat enonce dans [|[ permet d'etablir des inegalites de stabilite Lipschitziennes pour une classe de 
problemes inverses. II permet notamment de retrouver les resultats de stabilite developpes dans 12 1 et [l| sans avoir 



recours a des arguments de quantification de resultats de continuation unique. L'auteur precise que Ton peut trouver 
l'idee originale developpee dans Q dans le cas particulier de la detection d'un obstacle se deplagant dans un fiuide a 
partir de mesures disponibles sur le bord du domaine. De plus, des theoremes abstraits du meme type mais avec des 
hypotheses differentes peuvent etre trouves dans fl3j ]. 



2 Introduction 

Let T > 0, ft C R d , with d 6 N*, be a Lipschitz bounded connected open set such that dCl — L; U Lo U T out and 

N 

r otl t = r^. We are interested in the inverse problem of identifying the Robin coefficient q defined on some non 

i=l 

accessible part of the boundary r out from available data on Lo for (u,p) solution of the Stokes system (TTJ). Such kinds 
of systems naturally appear in the modeling of biological problems like, for example, blood flow in the cardiovascular 
system (see [ll[ and or airflow in the lungs (see Q). For an introduction on the modeling of the airflow in 
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the lungs and on different boundary conditions which may be prescribed, we refer to [8[. Similar inverse problems 

have already been studied in the stationary case in [f|, [j| and 0. In [B| and 0], a logarithmic stability estimate is 

obtained, whearas a Lipschitz stability estimate is established in [9( under the a priori assumption that the Robin 

coefficient is piecewise constant on T out . In each cases, the measurements involved in the stability estimates are the 

dp 

velocity u, the pressure p and the normal derivative of the pressure -— on F C Tq. The case of the nonstationary 

Stokes system has been addressed in [5[ in the particular case where the Robin coefficient does not depend time. The 
idea, introduced in [3j in the case of the Laplace equation, is to extend the stability estimate valid for the stationary 
problem to the nonstationary problem by using an inequality from the theory of analytic semigroups. This leads to 
infinite time measurements. 

The originality of the Lipschitz stability estimate presented in this Note is multiple: on the one hand, we 
obtain a stability estimate valid for the nonstationary Stokes system in finite time with a time-dependent Robin 
coefficient and secondly, the only measurement involved in the stability estimate is the velocity u on (0,T) x L, with 
r C To. In addition, the set of admissible Robin coefficients is more general than in Robin coefficients are not 
necessarily piecewise constant but belong to some compact and convex subset of a finite dimensional vectorial subspace 
of the set of continuous functions. Finally, we relax the regularity assumptions needed both on the boundary of the 
domain 57 and on the flux g. 

To be more precise, we introduce some notations. 

Notation 2.1. We denote by 

£+ ((0,T) x T out ) = {qe L°°((0,T) x T out ); 3m > 0, q > m a. e. on (0,T) x F otlt }, 

and 

C°(0,T;C° pc (T out )) = {q : (0, T) x T out -> R; ?| (0 ,T)xr 4 G C°((0,T) x F,) for 1 < i < N}. 

The main result of this Note is summarized in the following theorem. 

Theorem 2.2. Let M G N*. Let V M be a subspace of C°{0,T;C° c (T out )) spanned by some M linearly independent 
functions and Km be any compact and convex subset of Vm PI L3?((0,T) x F ont ). Let T C Tq be a nonempty open 
subset of the boundary, uq G be such that div uq = in fl and g € -ff^O, T; L 2 (Tq)) be such that g(t) is not 

identically zero for all t 6 (0,T). Let {ii]~,Pk) be the weak solutions of system ([T]) with q = <j% € Km for k = 1,2. 
Then, there exists a constant C > such that 

\Wl - <?2|U~((0,T)xr o „ t ) < C\\ui - M2||L 2 ((0,T)xr) 

The proof of Theorem 12.21 is based on an abstract theorem proved by L. Bourgeois in [6J that we recall in 
Theorem 13.41 and rely on the fact that the application which, to a Robin coefficient maps the velocity on (0, T) X T with 
r C Fo is injective, of class C 1 and its derivative is also injective. The result stated in Q establishs Lipschitz stability 
estimates for a class of inverse problems. For instance, it allows to find again the stability results developed in [12] 
and [l[ without resorting to quantification of unique continuation results. The author points out that one can find the 
original idea developed in 0] in the particular case of the detection of a moving obstacle in a fluid from measurements 
available on the boundary of the domain. Moreover, abstract theorems of the same type but with different assumptions 
can be found in The sequel of this paper is organized as follows. We present in Section [3] some preliminary results 
which will be useful to prove Theorem 12.21 Then the proof of Theorem 12.21 is given in Section @] 



3 Preliminary results 

In the section, we state results which will be useful in the proof of Theorcm l2.2l We begin by stating regularity result 
for a slightly more general Stokes system than system ([1} (we add non homogeneous Robin boundary condition on 
(0,T) x r oMt ): 



dtu — Am + Vp 


= o, 


in (0,T) x O, 


div u 


- o, 


in (Q,T) x n, 


u 


= o, 


in (0,T) x r,, 


d v u — pv 


= 9, 


on (0,T) x r , 


d v u — pv + qu 


= «j 


on (0,T) x r out 


u(0,.) 


= uq, 


in f2. 



(2) 
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Proposition 3.1. Assume thatg € H 1 (0,T;L 2 (r Q )), k G H 1 (0,T;L 2 (T out )), u G i2" x (fi) be such thatdiv u = in 
and? € L5°((0,T)xr oM t). T/ien, system (J2J) /ias a unique solution which belongs to L 2 (0,T; ff 1 (f2))ni/ 1 (0, T; L 2 (fi))x 
L 2 (0,T; L 2 (f2)). Moreover, there exists C > 0, independent of q, such that the following inequality holds 

IMlL2(0,T;ffi(^)) - C(ll u o||i3(n) + ll.9llL 2 (0 : T;L2(r )) + ll«llL 2 (0,T;L2(r out )))- 

Proof of Proposition \3.1\ The proof is mainly contained in the appendix of [f| . The main difference here is that we 
work with non homogeneous Robin boundary conditions which leads to slight modifications. □ 

Remark 3.2. Note that due to the mixed boundary conditions, the fact that dtu G L 2 (0, T; L 2 (Q)) does not imply that 
( U ,p)eL 2 (fl,T;iJ 2 (fi))xI 2 (0,T;ff 1 (l])) 

The following Proposition 13.31 concerns the identifiability of the inverse problem we are interested in. 

Proposition 3.3. Let T C To be a nonempty open subset of the boundary, uq € H (fl) be such that div uq — in fl. 
Assume that g 6 H 1 (0, T; i 2 (ro)) is such that g(t) is not identically zero for all t G (0,T). Let (uk,Pk) be the weak 
solutions of system (JTJ) with q = G C (0, T; C®jT out )) for k — 1,2. We assume that u\ = «2 on (0,r) x T. TTien 
9i = 92 ok (0,T) x r out . 

Proof of Proposition \3.S[ The proof is based on the unique continuation result for the Stokes system proved by C. Fabre 
andG. Lebeauin [ioj. Thanks to the previous proposition, (u,p) G L 2 (0, T; _ff 1 (il))n_ff 1 (0, T; L 2 (il)) x L 2 (0, T; L 2 (£7)), 
which is enough regularity to prove a similar result to Corollary 3.2 in [5]. Then, we proceed exactly as in the proof 
of Proposition 3.3 in [H[ where the proof is done in the particular case when the Robin coefficient does not depend on 
time, by arguing by contradiction. □ 

As announced previously, the proof of our Lipschitz stability estimate is based on a theorem proved by L. 
Bourgeois in Q which establishs Lipschitz stability estimate in an abstract framework for parameters defined on some 
finite dimensional subspace of the set of the continuous functions. For the sake of completeness, we state this theorem 
below. 

Theorem 3.4. Let (V, \\ \\v) and (H, \\ \\h) be two Banach spaces. Let U be an open subset of V and Vm a, finite 
dimensional subspace of V of dimension M . Let Km be a compact and convex subset of Vm H U . We consider a 
mapping T : U — > H which satisfies the following assumptions: 

1. T : Vm H U — > H is infective, 

2. T : U — > H is C 1 : T is differ entiable in the sense of Frechet at any point x G U, the Frechet derivative being 
denoted dT x :V^H and the mapping i£(/-> dT x G C(V, H) is continuous. 

3. For all x G Vm H U , the operator dT x : Vm — ?> H is injective. 

Then, there exists C > such that \/x,y G Km, \\x — y\\v < C||r(cc) — T(u)\\h- 

4 Proof of the main result 

In this section, we establishes the proof of Theorem l2.2l The proof consists of applying the abstract Theorem l3.4l with 
V = L°°((0,T) x r out ), H = L 2 ((0,T) xT), U — L9°((0,T) x T out ). We consider the operator 

T : L^((0,T)xT out ) L 2 ((0,T)xr) 
q -4 u\ r 

where (u,p) is solution of system (TT]). We are going to prove that: 
1. T\v M nu is injective, 
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2. T is differentiable at any point q £ L3°((0, T) x r out ) and its Frechet derivative is the operator 



where (vh,Th) is solution to 



dT q :L°°((0,T)xT o 



d t v h - Av h + Vr h = 0, 

div Vh = 0, 

Vh = 0, 

c^Wfi - T h v = 0, 

d v v h -T h v + qv h = —hu, 

«h(0,.) = 0, 



L 2 ((0,T) x r) 



in (0,T) x O, 
in (0,T) x Cl, 
in (0,T) x r,, 
on (0,T) x r , 
on (0,T)xr mt , 
in Cl, 



(3) 



where u is solution to system (JXJ) . Moreover, the mapping 



dT:i~((0,T) xT 



-»• r(L~((o, t) x r out ), £ 2 ((o, t) x r)) 



(4) 



the operator dT q : Vm -> £ 2 ((0,T) x T) is injective. 



is continuous. 
3. For all q £ V M n L£((0,T) x r otlf ) 

Step 1 is a direct consequence of Proposition I3~3l Let us prove step 2. Let q, h £ ((0,T) x F out ) and (u,_p) (resp. 
(uh,Ph)) be the weak solution of system ([1} associated to q (resp. to q = q+h). We denote by (wh, i^k) — (uh~ u )Ph—p) 
which is solution of the following Stokes system: 



d t w h - Aw h + V7r/j 


= o, 


in (0,T) x Cl, 


div Wh 


= o, 


in (0,T) x Cl, 


w h 


= o, 


in (0,T) x r,, 


dyW h -TThV 


= o, 


on (0,T) x r , 


dvWh -n h is + qw h 


= -hu h , 


on (0,T)xr mt 


w h (0,.) 


= o, 


in Cl. 



Let Mi > be such that ||ito||z 2 (n) + II<?IIl 2 (o,T;L 2 (To)) ^ Then, thanks to Proposition 13.11 there exists C > 
such that 

\\uu - u\\ L 2( T . H i^ < C\\h\\L°°((p,T)xr ou t)\\uh\\L a (o,T-,H' i -(n.)) < C(Mi)||/i||ioo(( 0jT ) x r oi4t ). (5) 

Now, let us consider (eh, Ph) — (uh — it — Vh,Ph —p — T~h), where (vh,Th) is solution to system First, we readily 
check that the operator h £ L°°((0,T) x T out ) — > Vh\r £ L 2 ((0,T) x T) is linear continuous. Secondly, (eh,Ph) solves 
the problem 



d t e h - Ae h + Vp h 


= o, 


in (0,T) x Cl, 


div eh 


= o, 


in (0,T) x Cl, 


e h 


= o, 


in (0,T) x r,, 


d v e h - PhV 


= o, 


on (0,T) x T , 


d„e h - PhV + qe h 


= -h(uh-u), 


on (0,T) x T out 


e h (0,.) 


= o, 


in f2, 



which implies, thanks to Proposition 13. II and inequality (JSJ), that 

||eh[U a ((o,r)i.ffi(n)) < C||^l|L-((o,T)xr out )ll u 'i _ u\\l^(o,T;H 1 (q)) < C(^i)ll /i lli~((o,T)xr out )' 

which proves that T is Frechet differentiable and dT q (h) — Vh\r- 

Let us prove now the continuity of the mapping dT defined in @. Let (vh,Th) (resp. (v h , r h )) be the solution 
of system §5§ associated to q (resp. q = q + l) and where (u,p) (resp. (u,p) — (ui,pi)) is the solution to system (JTJ) 
associated to q (resp. q = q + I) . We have that (v l h — Vh,p l h — Ph) is the solution to the following Stokes system: 



dt(v{ - v h ) - A(v l h - v h ) + V( P l h - Ph ) 
div (v l h - v h ) 



Vh 



dv(v l h -v h )- (p l h -p h )v 
d v (v{ - v h ) - (p h - Ph)v + q(v l h - v h ) 
(v l h -v h )(0,.) 



-lv l h - h(ui - it), 



in (0,T) x Cl, 
in (0,T) x Cl, 
in (0,T) x r,, 
on (0,T) x T , 
on (0,T)xr mf , 
in Cl. 
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This implies, thanks to Proposition 13. 11 

\\ v h ~ Vh\\L2(0,T;H l (n)) < C|UI|z,°-((0,T)xr o „ t ) H^i |U 2 (o,r;// 1 (")) + ( ^II^IU oo ((0,T)xr out )||M; - "|U 2 (0,T;Hi(O)), 

which leads to, applying again Proposition l3.il and inequality ([5]) with h = I: 

IK - u h ||x,2( 0iT; Hi(n)) < C(Mi)\\l\\ L o°(( 0tT ) x r out )\\h\\L°°((o,T)xr ou1 ,), 
where C is uniform with respect to h and I. Otherwise, we have proved that 

\\\dT q+l - dT q \\\ < C[[Z|| L oa ( (o,T)xr out ), 
where | | denotes the operator norm. Thus the mapping dT is continuous. 

It remains to prove Step 3. Let q G ((0, T) x r out ) n Vm ■ Assume that h s Vm is such that Wh|(o,T)xr = 0. 
Then, since (v/t|(o,T)xri {d u Vh — '?7i^)|(o.T)xr) — (0,0), we obtain from unique continuation result that Vh = in 
(0,T) x f2 and then hu = on (0,T) x r out . We conclude that h = by contradiction, exactly as for the injectivity 
of the mapping T (see [5|). 
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